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A fusion reactor requires plasma pre-heating before the rate of deuterium-tritium fusion reactions 
becomes significant. In ITER, radiofrequency (RF) heating of 3 He ions, additionally puffed into 
the plasma, is one of the main options considered for increasing bulk ion temperature during the 
ramp-up phase of the pulse. In this paper, we propose an alternative scenario for bulk ion heating 
with RF waves, which requires no extra 3 He puff and profits from the presence of intrinsic Beryllium 
impurities in the plasma. The discussed method to heat Be impurities in D-T plasmas is shown to 
provide an even larger fraction of fuel ion heating. 


I. INTRODUCTION 

In future fusion devices, radio frequency plasma 
heating with waves in the ion cyclotron range of 
frequencies (ICRF) is the only method capable of 
providing a significant fraction of bulk ion heating. This 
is due to the fact that electron cyclotron resonance 
heating, by definition, deposits wave energy to electrons. 
The same holds for the neutral beam injection systems 
because they will need to operate at rather high beam 
energies (~ 1 MeV) in order to reach the high-density 
plasma core. 

As discussed in [1-3], preferential bulk ion heating 
during the ramp-up phase has a number of advantages. 
These include an improved control over the path to the 
burn phase and keeping away from the region where the 
interaction between thermal ions and thermal electrons 
is weak. As a result, reaching the operational point 
with higher Q can be done faster with less auxiliary 
heating power than for an electron heating scenario [1], 
In addition, Tore Supra experiments demonstrated that 
higher levels of ion heating were accompanied by an 
improved energy confinement [2]. 

Applying ICRF power - somewhat counter-intuitively 
to its name - does not necessarily result in the dominant 
heating of bulk ions. In fact, depending on the chosen 
operational conditions (e.g., ICRF frequency, antenna 
phasing, plasma composition), the incoming RF power 
can be directly absorbed by various ion species, as well 
as by electrons. It is the latter condition one needs to 
fulfill in order to drive non-inductive current with ICRF 
system [4, 5]. However, channeling the RF power to be 
mostly absorbed by ions is not sufficient for bulk ion 
heating. In addition, one needs to avoid the formation 
of a high-energy ion tail due to ICRF. The critical 
energy F cr ; t of fast ions, at which the rate of loss of 
energy to the electrons and to the ions is equal, is 

given by F cr it = 14.8Af ast T e (J2i X l Zf/A l ) 2/ ' i [6]. Here, 
Af as t is the atomic mass of the energetic ions, T e is the 


electron temperature, = rii/n e , Zi and A are the 
concentration, the charge state and the atomic mass of 
the thermal ion species, respectively. In order to achieve 
the resultant dominant heating of bulk ions, the energies 
of RF-accelerated minority ions should stay in the range 
of Fcrit or below. 

ICRF scenarios relevant for the D-T plasmas were 
discussed in detail in Refs. [1, 3, 4, 7], with the emphasis 
to ITER in, e.g., [8, 9]. Many of these scenarios were 
studied experimentally during the past D-T experiments 
in JET and TFTR tokamaks [10-14], Currently, the 
second harmonic heating of tritium ions (w = 2uj c t) 
is considered as the main ICRF scenario for the ITER 
burn-phase plasmas, and the fundamental heating of 
a small fraction of 3 He ions (oj = w c 3 He ), additionally 
puffed into the plasma, for the ramp-up low-temperature 
heating stage of the pulse. In JET, a fourfold increase of 
the neutron emission was observed when A'[ 3 He] = 4% 
was injected [11]. The higher reactivity observed 
with 3 He minority heating was due to the higher Xj 
reached. However, currently the supply of 3 He reduces 
and the industrial demand of this gas is progressively 
increasing [15]. 

In this paper, we introduce and discuss a promising 
alternative bulk ion heating scenario for the ramp-up 
phase in D-T plasmas. We will mainly focus on the 
applicability of the proposed scheme for ITER. Instead 
of using additionally puffed 3 He ions, we suggest to tune 
the RF system to heat predominantly intrinsic Beryllium 
(Be) impurities, which will be naturally present in ITER 
plasmas. Because of the larger atomic mass of Be, such 
ICRF heating is shown to provide an even larger fraction 
of bulk ion heating. Finally, we discuss how the current 
design of the ICRF system in ITER can be adapted to 
exploit the advantages of beryllium ICRF heating to their 
maximum. The proposed method is not restricted to 
the use of °Be only, but other impurities with a similar 
charge-to-mass ratio, such as 'Li, 22 Ne, 40 Ar, etc. can 
be used as well. This is discussed in the last section of 
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the paper and is illustrated with a few examples. 


II. THREE-ION (Be)-D-T ICRF SCENARIO 


ICRF heating relies on launching the fast 
magnetosonic wave (FW) into the plasma by external 
antennas, typically located at the low field side (LFS) 
edge of the plasma. The FW propagation across the 
plasma is fairly well described by the well-known 
dispersion relation [16] 


(<T - nfMen - n\) 
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where nju = c/cj_,||/ w is the perpendicular/parallel FW 
refractive index, and the tensor components es, and 
cr, are those given by Stix [17]. 

The FW power is absorbed by ions, when the wave 
crosses the ion cyclotron (IC) resonance (ui = oj C i ) and 
IC harmonic (w = Nui a , N > 2) layers. Electron 
absorption is also possible via a combination of electron 
Landau damping and transit time magnetic pumping, 
especially at high electron beta. Electron heating can 
also occur via mode conversion of the incoming FW 
to shorter wavelength modes, which takes place at the 
ion-ion hybrid resonance(s) in multi-ion plasmas [18]. 

In general, the FW is elliptically polarized and its 
electric field can be decomposed as a sum of two 
components: the left-hand polarized component E+, 
which rotates in the sense of ions in the magnetic field, 
and the right-hand polarized component that is 
aligned with electron rotation. In his seminal paper [19], 
Stix showed that the absorbed RF power by thermal ions 
is almost merely due to the left-hand E + component. 
One should note here that it is the plasma, rather than 
the ICRF system, which defines the wave polarization 
and imposes the ratio between E + and E_ [19] 
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Though the characteristic ny, which appears in Eqs. (1) 
and (2), can be varied by changing the ICRF antenna 
phasing and operational frequency, the main contribution 
to Eq. (2) generally comes from the tensor elements. 

One option for ion absorption in fusion plasmas is 
applying IC harmonic heating (TV > 2), \E + /E_\ ~ 
(.N — 1)/(N + 1). This damping mechanism is, however, 
a finite Larmor radius effect and hence not appropriate 
for the beginning of the heating phase in a reactor 
(low densities and low temperatures). Another option 
known as minority heating in two-ion plasmas [19] 
is to heat a small fraction of resonant minority ions at 
their fundamental IC resonance. This scenario typically 
shows the best performance at minority concentrations 
of A m i no = rimino/n e ~ 3 — 10% and has been routinely 
used for plasma heating. Though minority heating has 
a quite strong single-pass absorption (SPA), especially 
for hydrogen minority heating in deuterium majority 
plasmas, the ratio of the left- to the right-hand 


polarization at the IC resonance of minority ions is 
limited [11] 
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and the FW is mostly right-hand polarized. Here, 
we introduced the notation Zi = (Z/A)i, which 

stands for the ratio of the charge state to the atomic 
mass for the ion species. By way of example, for 
(H)-D heating \E + /E_\ « 1/3, which is enough to 
obtain an efficient single-pass absorption at Ah ~ 
5% (typical H concentration used in the experiments). 
Yet, two-ion ICRF minority heating at much lower 
minority concentrations (e.g., A m i no <C 1%) is usually 
characterized by very poor wave absorption and for this 
reason is not often used in practice. 

In our recent paper [20], we presented a method 
to maximize the E + component and, subsequently, to 
achieve a strong ion absorption at very low minority 
concentrations (~ 1% or lower). This enhanced ion 
absorption requires the presence of three ion species in 
a plasma (with a different Z/A ): two main ion species 
(X and Y) and a third resonant ion species (Z), present 
at a very small concentration. The key idea is to 
locate the left-hand polarized fast wave L-cutoff, which 
is defined by the condition cl = and which largely 
enhances the associated RF field component E +, close to 
the fundamental cyclotron resonance of the minority ion 
species Z. This matching can be achieved by adjusting 
the density ratio between the two main ion species Y:X. 

Note that not all ICRF scenarios with three ion species 
can lead to E + enhancement. One of the necessary 
conditions is that 

mm{Z 1 ,Z 2 } < Z 3 < max{Zi, Z 2 }, (4) 

i.e. the resonant minority ions, one aims to heat, 
should have a cyclotron resonance located between the IC 
resonances of the two bulk ion species. The main ICRF 
scenario for bulk ion heating in D-T ITER plasmas, viz. 
( 3 He)-DT, does not belong to this list since Z 3 > Z\ )2 . 

Tungsten and beryllium have been chosen as the 
plasma-facing components for ITER. Accordingly, the 
carbon wall was replaced with the ITER-like wall in 
JET since 2009. Therefore, JET and ITER plasmas 
will unavoidably include a non-negligible amount of Be 
impurities. For example, Ae e = 1 — 2% is estimated as 
the background Be level in JET plasmas [21, 22]: the 
measured core Z e g « 1.2 yields Ase — 1.4% for the 
residual carbon content of about Ac ~ 0.1% (note that 
Z e g « 1.1 yields Aee ~ 0.6%). A charge-to-mass ratio of 
the fully-ionized beryllium ions {A = 9, Z = 4) satisfies 
the condition of Eq. (4), Zt < Zs e < -Zb, and therefore 
(Be)-D-T scenario belongs to the set of three-ion ICRF 
heating scenarios. Though in Ref. [23] it was correctly 
outlined that the application of D-T mode conversion 
heating could be hampered due to the Be presence, our 
recent results [20] show that the impact of the three-ion 
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FIG. 1. (a) Location of the IC resonances in the ITER plasma at the full magnetic field (Bo = 5.3 T). The official ITER 

frequency range for the ICRF system / = 40 — 55 MHz is highlighted, (b) The normalized L-cutoff frequency, <I>l = wl/wcH as 
a function of tritium and beryllium concentrations in a (Be)-D-T plasma. 


ICRF scenarios can be reverted and such scenarios have 
a great potential for fusion research (fast-ion generation 
and bulk plasma heating). Below, we explore the features 
of the (Be)-D-T ICRF scenario in ITER-like plasmas and 
compare the potential of bulk ion heating for this scenario 
with the more traditional ( 3 He)-DT scheme. 

III. (Be)-D-T HEATING IN ITER 

Figure 1(a) shows the location of the main IC 
resonances in ITER (Rq = 6.2 m, a = 2.0 m) as a function 
of the ICRF frequency for the nominal magnetic field 
Bo = 5.3 T. In ITER, two ICRF antennas mounted 
in mid-plane ports are envisaged to operate within the 
range / = 40 — 55 MHz [24, 25]. Such a frequency 
range for ITER provides central N = 1 and N = 2 
heating for 3 He and T ions at / « 53 MHz, and 
fundamental IC heating of D ions at / ~ 40 MHz (for 
(D)-T pulses). We also depict the location of the IC 
resonance of Be impurities (red solid line) in Fig. 1. One 
can see that the distance between the Be and D cyclotron 
resonances is about 70 cm, and a frequency only slightly 
below the official minimum ICRF frequency for ITER 
/ « 36 — 37 MHz is required to locate the Be resonance 
close to the plasma center. A frequency / = 37 MHz 
is used throughout the paper, unless otherwise stated. 
According to the arguments presented in the last section 
of this paper, the ICRF system in ITER can operate 
at this frequency, without any dramatic changes to the 
present design. As will also be discussed below, operation 
at / = 40 MHz can still provide dominant Be heating at 
the price of yielding reasonably off-axis power deposition. 

Depending on the impurity concentration, one can 
distinguish two different regimes of the three-ion ICRF 
heating. If the concentration of Be ions is extremely 
small (XB e <C 1%), the optimal T concentration for 
maximizing Be absorption can be found from Eq. (7) of 
Ref. [20], viz. X| » ( Z Be - Z T )/(Z B - Z T ) » 67% 


(D:T ~ 1:2). At such Xt, the L-cutoff in D:T plasmas 
(see dashed-dotted line in Fig. 1(b)) and the cyclotron 
resonance of Be impurities intersect. However, in 
three-ion Be-D-T plasmas with a finite X Be: the L-cutoff 
can not cross the IC resonance of Be. This is clear 
from the equation for the normalized L-cutoff frequency 
(w L = wl/w c h), which can be written as follows [26] 


Xt Xd 4Xg e 

WL — Zt Wl — -Zb Wl — Z Be 


(5) 


Note that any solution of Eq. (5) can be translated 
into the corresponding radial coordinate of the L-cutoff 
via i?L ~ R 0 wl 15.25 Bo(T)//(MHz). In fact, as follows 
from Fig. 1(b), in the vicinity of u> = w c , Be, the curve for 
wl splits in two, and generally there are two meaningful 
solutions of Eq. (5) at a given Xt and X Be . Far from the 
region of the resonant crossing (Xt ~ X.f), they can be 
approximated with expressions 


WLl ~ z Be + 


4 (Z D - Z Be ) X.p X B e 


X4,-X t 


u> L2 ~ -Zt + (Z B — Zt)Xt — 


4(-Zp — -Z’Be)XT Xse 
X* -X T 


( 6 ) 


In the vicinity of Xt ~ X^, these approximations 
are not valid. For very low Xs e -C 1%, the gap 
between two solutions wli and (I>L 2 is small, and the 
optimal conditions for the ICRF heating occur at the 
tritium concentration X^., at which the two asymptotics 
intersect. We have depicted such conditions as ‘Regime I’ 
in Fig. 1(b), and this heating regime has been suggested 
for fast-ion generation with ICRF and studied in detail 
in [20], 

However, for higher Be concentrations (XB e ~ 1%), 
at Xt — X£ the gap between wli and w Cj Be is already 
quite large. As a result, the distance between the 
region with the enhanced E + and Be cyclotron resonance 
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FIG. 2. Single-pass absorption coefficients as a function of the ICRF frequency for the ITER, ramp-up phase as computed by 
TOMCAT (B 0 = 5.3 T, D:T=1:1, T 0 = 4 keV, n e0 = 6 x 10 19 m -3 , A;jj ant) = 3 m _1 ): (a) ( 3 He)-DT scenario (A 3H e = 2%, 
Ase = 1%); (b) (Be)-D-T scenario (no A3 He , Ase = 1%). 


exceeds the Doppler width of the IC resonance, and 
minority ion absorption at such conditions significantly 
decreases. Thus, as already noted in [20], the optimal 
tritium concentration for maximizing Be absorption 
shifts towards lower values of At- These conditions, 
which are schematically depicted as ‘Regime II’ in 
Fig. 1(b), correspond to the conditions of (Be)-D-T 
heating in ITER and will be studied in this paper. 

We start with a comparison between the ( 3 He)-DT and 
(Be)-D-T ICRF scenarios by evaluating the single-pass 
absorption coefficients (a fraction of the incoming RF 
power absorbed by plasma species as the wave propagates 
from the LFS edge to the high field side edge) as 
a function of ICRF frequency. These are computed 
with the TOMCAT code [27] and are depicted in 
Fig. 2. For both scenarios, the background X Be = 1% 
and an optimal D:T=1:1 ratio is assumed. For 
the 3 He minority scenario, shown in Fig. 2(a), an 
additional A3 He = 2% is considered. Because of the 
plasma dilution, the concentration of fuel D-T ions for 
this scenario is somewhat lower than for Be minority 
heating (Ad,t = 46% vs. Ad,t = 48%). We focus on 
the ITER ramp-up phase, when the plasma density 
and temperatures are rather low and the rate of D-T 
fusion reactions is small, such that the contribution 
of alpha-particles to the power balance is negligible. 
As baseline conditions, we consider To = 4 keV, 
n e o = 6x 10 19 m -3 , ntor = 27 (fcjj ant ' ) = 3 m _1 ). Such 
a value of the FW parallel wavenumber is representative 
for the [0; 0; n ; 7r] toroidal phasing of the ITER ICRF 
antenna [8, 28]. An influence of the antenna toroidal 
spectrum on ICRF heating scenarios in ITER was 
studied in [9, 28]. For the ( 3 He)-DT scenario, two 
modelling approaches (full-spectrum computations and 
retaining only the dominant toroidal mode number) were 
shown to give almost identical results (see Table II of 
Ref. [28]). Hence, for proof-of-principle computations 
which we present in this paper, we have adopted a simpler 
approach and considered only the dominant wavenumber. 


As follows from Fig. 2, for the considered conditions 
the single-pass absorption for the 3 He minority scenario is 
P 3 He ~ 40% at / « 50 MHz. The single-pass absorption 
for the Be heating scenario can be stronger p Be « 90%, 
if operating at a lower frequency / ~ 37 MHz. Note 
that p 3He ~ 40%, which is computed for the ( 3 He)-DT 
scenario, does not mean poor ICRF heating. Typically, 
such pspa is still high enough for having good multi-pass 
ICRF performance (pspa = 40—50% reflects that the FW 
needs to pass the resonant layer several times before the 
power is absorbed). Note that operation at / ~ 43 MHz 
(secondary peak of Be absorption in Fig. 2(b)) is not 
very promising since then the RF power is deposited 
appreciably off-axis ((r/a )Be ~ 0.5, cf. Fig. 1(a)). 

There is another criterion to consider: one needs 
to make sure that the ICRF scheme in question is 
robust with respect to changes in the D:T ratio and Be 
concentration. Since the results shown in Fig. 2 were 
computed with a ID code, which neglects a number of 
effects intrinsic for the tokamak geometry such as a finite 
poloidal magnetic field, a more sophisticated modelling 
with the 2D full-wave codes TORIC [29, 30] and EVE [31] 
was done to increase the degree of realism. Note that 
contrary to the ID computations, in 2D modelling (since 
the tokamak vessel acts as a Faraday cage) all the 
RF power is eventually absorbed in the plasma, and 
the computed absorption coefficients reflect the relative 
strength of various damping mechanisms. For the 2D 
computations we adopted a Shafranov shift of 10 cm 
and applied Bq(Ro = 6.2 m) = 5.4 T such that the 
magnetic field at the magnetic axis is 5.3 T as for the 
ID computations. There is generally a good agreement 
between TORIC and EVE results. Figure 3(a) shows 
that damping of the RF power on Be impurities is the 
dominant absorption channel (pee > 60%) for a wide 
range of tritium concentrations, namely At = 20 — 60%. 
This is in line with the arguments behind Fig. 1(b) and 
discussions above. The remaining incoming RF power is 
almost equally split between fuel D ions and electrons. 
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FIG. 3. Beryllium impurity ICRF heating scenario, (Be)-D-T: (a) Absorption coefficients as a function of tritium concentration 
(Xae = 1%), as computed by TORIC, EVE and TOMCAT; (b) Radial dependence of the real part of fc^ FW (A'ue = 1%, 
.Yd = Xt = 48%); (c) p a bs vs. beryllium concentration (D:T=1:1). 


One can also notice a modest oscillatory behavior 
of pBe in Fig. 3(a). TORIC computations predict 
that the maxima of Be absorption are reached at 
Xt ss 22%, 34%, 44% and 58%, and minima at Xt ~ 
28%, 39% and 54%, respectively. This occurs due to 
a constructive/destructive interference effect. Similar 
to the structure considered in [32-35], for (Be)-D-T 
plasmas there are two MC layers located in the plasma 
(see Fig. 3(b)). As noted, Eq. (5) has two meaningful 
solutions: the first L-cutoff Rbi is located close the IC 
resonance of Be ions, R Be (to the lower magnetic field 
side of that) and the second R B 2 to the higher held 
side of I?Be- hr contrast to ‘Regime I’ heating, when 
the impurity concentration is too small to form the two 
well-separated MC layers, for this regime the RF power 
is not very efficiently absorbed, when the FW passes 
through the first MC layer and Be cyclotron resonance 
(pi ~ 30%). The transmitted FW (T — 40%) reaches the 
second MC layer and is almost entirely reflected back at 
I?L 2 - As a result, the interference of the two reflected 
waves determines the resultant reflection and absorption 
coefficients. Analytical formulae describing the 
constructive/destructive interference effect for minority 
heating regime were derived in [36]. According to 
Eqs. (5) and (6) of [36], the minority ion absorption varies 
between pi.min ~ 13% and pi, ma x ~ 99%, depending 
on the phase difference A(j> between the two reflected 
waves. An average single-pass absorption coefficient is 
then pi = pi + T(1 — T) ~ 54%. 

The two MC layers shown in Fig. 3(b) have essentially 
different radial width. Whereas for the first layer 
Amci — 1 — 1-5 cm, the second MC layer is much wider 
AmC 2 — 10 cm and is, thus, non-transparent for the wave 
propagation. As the concentration of T ions is gradually 
decreased from XJ, the distance between I?l 2 and the 
first MC layer increases (see Eq. (6)) 

A R 12 ~ [(Z Be - Z T ) - (Z D - Z T )X T ] RBe/Z Be « 

« (1 — 3X t / 2) l? 0 /4. 


Every new maximum/minimum shown in Fig. 3(a) 
marks a case, when the phase term $ ~ ^±.fwAI?i 2 
has changed by d4> = 7 t/ 2, i.e. the radial distance R 12 
has increased roughly by a quarter of the average FW 
perpendicular wavelength. Naturally, such oscillations 
are much more pronounced in ID (TOMCAT) 
computations, and the modulation amplitude for 
2D p a bs-computations is much smaller. Since 
k±.FW - {u p R/c)Z Be (J2i X i A i) 1/2 ~ 0.7(w pH /c) with 

w p h = (47m e e 2 /mH) 1//2 , a characteristic tritium 
variation, marking the transition from the maximum of 
p Be to the minimum and vice versa, has been derived 

SX T ^-p- «0.27/(VS^i2o(m)). (7) 

o/C_[_jT10 

Here, n e ,20 is the plasma density expressed in the units 
IO 20 m -3 and Rq is the major radius of the machine in 
meters. 

Accurate measurements of the concentration of 
Be impurities in the plasma center are not readily 
available. Though we fixed X Be = 1% in our previous 
computations, one should consider this value only as a 
rough estimate. Figure 3(c) depicts the sensitivity of the 
2D absorption coefficients as a function of Xse- The 
concentrations of fuel D and T ions for every numerical 
run were adjusted accordingly to obtain an optimal 
D:T=1:1 mixture. Our results show that the RF power 
is absorbed mostly by Be impurities within a fairly wide 
range XB e ~ 0.1% — 2%. For the considered conditions, 
the Be absorption has peaks around XBe ~ 0.5% and 
1.2%, and tends to decrease at higher concentrations. 
At XBe > 2.4%, electron absorption starts to be equally 
important. Should ITER plasmas include high Be 
levels XBe > 2.5 — 3%, this will result in a reduced 
Be absorption. However, under these conditions at least 
20% of the output fusion power will be lost due to the fuel 
dilution. A general conclusion, which one can draw from 
studying Figs. 3(a) and (c), is that for the given ICRF 
frequency / = 37 MHz Be absorption is the dominant 
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FIG. 4. An effect of extra impurity species with Z/A « 0.43 — 0.45: absorption coefficients as a function of (a) 7 Li and (b) 22 Ne 
concentrations. Lines with symbols correspond to the case Xs e = 1% (data for Ane = 2% case is shown with lines only). 


TABLE I. Location of IC resonant layers u> = to c i for bulk 
ion and impurity species, xic = Ric — Ro ( Ro = 6.2 m, 
a = 2.0 m). 



Ion species 

T 

'Li 

Be 

22 Ne 

D, 2U Ne 

37 MHz 

xic (m) 
xio/a 

-1.69 

-0.85 

-0.40 

-0.20 

-0.19 

-0.09 

-0.05 

-0.03 

0.57 

0.28 

40 MHz 

xic (m) 
xic/a 

-2.03 

-1.01 

-0.84 

-0.42 

-0.64 

-0.32 

-0.51 

-0.26 

0.06 

0.03 


channel of the RF absorption for a wide range of realistic 
plasma conditions in ITER. 

The dashed-dotted lines shown in Fig. 3(a) and (c) 
depict the absorption coefficient by Be impurities 
computed with TORIC for / = 40 MHz. It is clear 
that Be absorption dominates over the other damping 
mechanisms for this ICRF frequency as well. 

The presence of high -Z impurities like tungsten has 
a small impact on the (Be)-D-T scenario. A typical 
charge state for W in the plasma core is Z ~ 40 — 50 [37] , 
and such impurities have (ZjA)^ « 0.2 — 0.3, which is 
very different to {Z/A) Be- On the contrary, impurities 
with a charge-to-mass ratio close to that of Be ions, e.g. 
7 Li , could have a non-negligible effect on RF power 
absorption. Figure 4(a) shows absorption coefficients by 
Be and 7 Li impurities computed by TORIC for D:T=1:1 
plasmas. As follows from this figure, the total absorption 
of RF power by two different impurity species pi mp = 
PBe+P 7 hi remains nearly constant. Since Z/A for Be and 
' Li ions is not identical, the IC resonant layers of these 
impurities are separated by a distance about 20 cm, and 
therefore the radial position of the maximum of impurity 
absorption is somewhat different (peaking at r/a — 0.1 
and 0.2, respectively). 

Another impurity species relevant for the studied 
ICRF scenario is the second most abundant isotope 
of neon, 22 Ne 10+ (natural abundance 9.3%). In fact, 
for / = 37 MHz the IC resonance of 22 Ne impurities is 
located closer to the plasma center than the Be resonance 
(see Table I). As a result, FW excited by a LFS antenna 
approaches the resonance of 22 Ne ions first and hence 
one can tune most of RF power to be absorbed by 


TABLE II. RF power split (direct heating and collisional 
redistribution) for the ( 3 He)-DT and (Be)-D-T ICRF 
scenarios. 



Electrons 

3 He/Be ions 

Bulk ions 

a) ( a He)-DT scenario: 

Direct RF heating 

9% 

76% 

14% 

Collisional redistribution (10 MW) 

29% 

5% 

62% 

Collisional redistribution (20 MW) 

43% 

4% 

50% 

b) (Be)-D-T scenario: 

Direct RF heating 

6% 

82% 

12% 

Collisional redistribution (10 MW) 

10% 

9% 

81% 

Collisional redistribution (20 MW) 

17% 

7% 

76% 


neon ions. This is clearly illustrated in Fig. 4(b), where 
a scan in A 22 Ne is made for two beryllium background 
levels. E.g., for Xee = 1%, adding as little as 0.02% 
of 22 Ne ions allows changing the dominant impurity 
absorption channel and bringing the absorption region 
even closer to the plasma center. Another stable isotope 
of neon, 21 Ne, is very rare (abundance 0.27%), and for 
the considered conditions not absorbing much of RF 
power. A preliminary conclusion one can draw from 
Figs. 4(a) and (b) is that having extra impurity species 
with a similar Z/A « 0.43 — 0.45 and IC resonant layers 

located close to Be resonance ( /a ~ 0.1, 

see Table I) is not hindering the here considered heating 
scenario. 

Comparison of the collisional power 
redistribution 

Finally, we compare the collisional power 
redistribution for the 3 He and Be minority heating 
scenarios. As outlined in the introduction, the critical 
energy, at which an equal amount of absorbed RF 
power collisionally goes to bulk ions and electrons, 
increases with the atomic mass of the energetic particles: 
E C rit [ 3 He] « 25T e and E cr ; t [Be] ~ 74T e , respectively. 
Hence, at comparable ICRF power densities, a stronger 
bulk ion heating is expected for the (Be)-D-T scenario. 

In terms of the fraction of RF power absorbed directly 
by various plasma species, the two ICRF scenarios are 
very similar. As follows from Table II, minority ions 
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( 3 He and Be, respectively) absorb about 75 — 80% of the 
launched power, < 10% of the RF power is deposited on 
electrons, and the rest 10 — 15% is absorbed by bulk D 
and T ions. Using the RF power densities evaluated by 
TORIC, the Fokker-Planck solver SSFPQL [38] was used 
to estimate the collisional power redistribution between 
the species. For Picrf = 10 MW of coupled RF power, 
the fraction of bulk ion D-T heating is 62% for the 
( 3 He)-DT scenario. This reduces to 50% for 20 MW 
of coupled RF power. The electron heating fraction is 
computed to be 29% and 43%, respectively. Note also 
that a finite amount of RF power eventually goes to heat 
minority 3 He and Be impurities. The computed results 
are quite similar to those reported in Ref. [9] (however, 
a higher plasma density and temperatures were used in 
those simulations). 

For the (Be)-D-T scenario, the fraction of bulk ion 
heating is larger. For 10 MW of coupled RF power, 
about 80% eventually ends up in bulk D-T ions and 
only 10% in electrons. Increasing the coupled RF power 
does not result in any significant reduction of pi. For 
-Picrf = 20 MW, fuel D and T ions still absorb 76% of 
the coupled power and the fraction of electron heating 
increases only up to 17%. 

Note that if extra natural neon (with natural 
abundances of isotopes) is puffed into the plasma, this 
can result in a somewhat higher electron heating fraction. 
For example, adding 0.5% of natural neon gas (A Z e g « 
0.5) brings only ~ 0.05% of 22 Ne ions. As shown in 
Fig. 4(b), 22 Ne ions can absorb a lot of RF power 
even at such very low concentrations. Subsequently, 
a tail of high-energy ions is expected to develop in 
the 22 Ne distribution function which can collisionally 
transfer a significant fraction of power to the electrons. 
Such a more complicated multi-impurity option deserves 
a further detailed study (this scenario is also numerically 
more challenging because of the multiple ion-ion hybrid 
layers present in the plasma) and the corresponding 
results will be reported elsewhere. 

IV. DISCUSSION AND CONCLUSIONS 

A big advantage of the ( 3 He)-DT heating scenario is 
that it gradually converts to w = 2 w c t heating, without 
any change in the ICRF frequency. Thus, 3 He is needed 
only during the ramp-up phase of the pulse and its puff 
can be switched off once the plasma temperature is high 
enough for second harmonic tritium heating to become 
relatively strong. Good performance of this scenario has 
been proven experimentally and if the availability of 3 He 
is not problematic, this scenario should be considered as 
the main option for ITER. 

The results of our studies show that there is a good 
backup option, especially for the very early heating stage 
of the pulse. The design of the ITER ICRF system 
currently envisages the use of two antennas [24, 25]. 
With that, one could imagine several ICRF strategies for 
the ramp-up phase. For example, an interesting option 
would be to start a pulse with one antenna operating 


at 37 38 MHz (or 40 MHz) and selectively heat Be 
impurities. When the plasma pre-heating with the first 
antenna gets efficient, one can switch on the second 
antenna at 53 MHz and start 3 He puffing to increase 
Ti further. The first antenna is then switched off and 
the ICRF system is configured to launch the RF power 
at a higher frequency later on (e.g., either for w = 2 w c t 
heating or for current drive to extend the pulse length). 
In any case, by using the first antenna at 37-40 MHz 
when Tj is low, one can reduce the total consumption of 
3 He during the ramp-up phase. 

For ITER operating at full magnetic field 
Bq = 5.3 T, the proposed (Be)-D-T ICRF scenario 
requires / ss 37 MHz to achieve heating close to the 
plasma center. At first glance, this frequency is 
not ITER-relevant since the official ICRF range is 
/ = 40 - 55 MHz [24, 25]. However, ITER RF 
generators are foreseen to deliver RF power in the range 
/ = 35 — 65 MHz [39]. Furthermore, according to the 
computed ICRF antenna performance in ITER, the 
reduction of the plasma coupling if operating at a lower 
frequency is reasonably small. If compared to the results 
for the officially approved frequency 40 MHz, the coupled 
power is lower by about 15% only for / = 38 MHz and 
by 20% for 37 MHz [24, 40, 41]. 

Yet, if the extension of the frequency range for ITER 
is not possible, operating at / = 40 MHz sets no 
major physics limitations for the proposed (Be)-D-T 
scenario. A drawback of such operation is the non-central 
power deposition, with the heating maximum shifted 
about 60 cm to the high field side ( r/a « 0.3 — 0.35). 
Also for this configuration, one can expect a stronger 
effect of parasitic alpha-particle absorption since at 
this frequency the IC resonance of D and 4 He ions is 
located centrally. Direct absorption of ICRF power by 
energetic alpha particles is a potential showstopper for 
the (Be)-D-T heating scenario at higher Ti. Using the 
recently upgraded version of the TORIC code [42], we 
have estimated that this absorption is reasonably small 
up to To ~ 10 — 15 keV: at To = 10/15 keV the fraction 
of alpha-particles in the plasma core is computed to be 
about 0.2%/0.5% and they absorb p a ~ 7%/19% (pse = 
72%/52%) of the incoming RF power, respectively. 
However, channeling of fusion alpha-particle power using 
minority ion catalysis [43] could potentially facilitate the 
application of Be impurity heating in D-T plasmas at 
higher T 0 . 

We want to outline a few examples of the three-ion 
ICRF heating in D-T plasmas already reported in the 
literature: 

(a) During the first D-T experiments in TFTR, reaching 
localized electron heating and driving non-inductive 
current via mode conversion, was found to be 
unsuccessful. The reason for that was attributed to 
the presence of lithium impurities in the plasma [23]. 
Numerical computations showed that most of the RF 
power was absorbed by 7 Li ions (A7 Li « 0.5%). This 
absorption was considered as parasitic, and extensive wall 


conditioning with 6 Li-enriched pellets was used later to 
eliminate this effect and recover normal MC heating. 

(b) For the JET pulse 42769, Start et al. reported 

~ 1.5%. It was estimated that the Be impurities 
were absorbing about 40% of the RF power [11]. Such 
a non-negligible Be heating was observed even though 
the IC resonance of Be impurities was located off-axis 
and although the used D:T ratio (Xd = 18%) was quite 
different to the optimal values we have computed. 

(c) In Ref. [3] (Fig. 7), dominant absorption of RF power 
by Be ions in D:T=1:1 ITER plasmas was computed for 
XB e = 1%, / = 40 MHz and = 4 m -1 . 

(d) In recent papers, where the ICRF scenarios for 
the activated phase of ITER were studied, a significant 
impurity absorption, mostly by argon, was reported 
[8, 9]. Because the ICRF frequency was chosen to locate 
the 1 2 3 4 5 6 7 8 He resonance centrally, the argon resonance was 
placed at the very high field side edge. Note that fully 
ionized argon ions have a charge-to-mass ratio (Z/A = 
0.45) very similar to that of Be ions. Again this impurity 
absorption was considered as a parasitic effect; however, 
in [9] it was outlined that this effect deserves further 
investigation with respect to the ICRF operation in 
ITER. 

As follows from the examples shown above and using 
the reported results, one can conclude that Be can be 
a very important species for future JET and ITER 
D-T experiments involving ICRF heating. In view of 
ITER, this scheme can be tested during the forthcoming 


DTE2 campaign at JET [44]. For JET operating at 
Bo = 3.6 T, ICRF heating of Be impurities in D-T 
plasmas requires the lowest frequency available for the 
A2 antennas / = 25 MHz. 

Finally, the proposed N = 1 impurity heating scenario 
in D-T plasmas is relevant for DEMO and future fusion 
reactors. In fact, there are several choices for the 
resonant absorber species in D-T plasmas. If for any 
reason Be can not be used as a wall material for the 
tokamak-reactor DEMO, this can be equally replaced 
with the other impurities with a similar charge-to-mass 
ratio, viz. ' Li 3+ (ionization energy for the Li 2+ Li 3+ 
transition is £i on . = 122 eV [45]), 22 Ne 10+ (£i on . = 

1.4 keV), 40 Ar 17+ (£ ion . = 4.1 keV) and 40 Ar 18+ (£ ion . = 

4.4 keV), etc. Note that the abundance of lithium is 
no problem for a fusion reactor, and argon and neon 
are likely to be available for gas puffing as these noble 
gases can be used for impurity seeding and controlling 
the deposition of the heat load in a machine. 

Summarizing the arguments, the reported results show 
that heavy intrinsic impurities with 1/3 < Z/A < 1/2 
can be an efficient absorber of the ICRF power in D-T 
plasmas, and such a resonant impurity RF heating can 
be used for increasing the bulk ion temperature T) during 
the ramp-up phase of the plasma pulse. 
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